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Preface 


The  evaluation  and  comparison  study  reported  herein  was  performed  by  the 
Hydraulics  Laboratory  (HL)  of  the  U.S.  Army  Engineer  Waterways  Experi¬ 
ment  Station  (WES)  during  July  through  September  1994.  The  U.S.  Army 
Engineer  District,  Vicksburg,  requested  WES  to  perform  this  work  and  to 
document  all  findings.  The  scope  of  this  project  consisted  of  comparing  an 
Acoustic  Doppler  Current  Profiler  (ADCP)  to  standard  discharge  measurement 
techniques  currently  in  practice  on  the  Mississippi  River.  A  comparison  plan 
was  developed  to  demonstrate  both  techniques  for  two  flow  conditions  on  the 
Mississippi  River  at  five  historically  monitored  ranges.  WES  collected  the 
ADCP  data,  while  a  private  contractor  under  contract  with  the  District  took  the 
standard  discharge  data.  Both  sets  of  data  and  the  comparison  results  are 
contained  in  this  report. 

The  study  was  conducted  under  the  direction  of  Messrs.  Frank  A. 

Herrmann,  Jr.,  Director,  HL;  Richard  A.  Sager,  Assistant  Director,  HL;  and 
William  H.  McAnally,  Jr.,  Chief,  Waterways  and  Estuaries  Division,  HL. 
Project  engineer  for  the  study  was  Mr.  Thad  C.  Pratt,  Prototype  and  Field 
Studies  Group,  HL.  Other  personnel  participating  in  the  data  collection  were 
Messrs.  Howard  A.  Benson  and  Christopher  J.  Callegan,  also  of  the  Prototype 
and  Field  Studies  Group.  Data  reduction  was  performed  by  Mr.  Pratt.  This 
report  was  prepared  by  Mr.  Pratt  with  assistance  from  Ms.  Clara  J.  Coleman, 
Technical  Support  Group,  HL. 

Dr.  Robert  W.  Whalin  was  Director  of  WES  during  the  publication  of  this 
report.  COL  Bruce  K.  Howard,  EN,  was  Commander. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication, 
or  promotional  purposes.  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 


Conversion  Factors, 
Non-SI  to  SI  Units  of 
Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units 
as  follows: 


1  Introduction 


The  Vicksburg  District  monitors  multiple  gaging  stations  along  the 
Mississippi  River  using  conventional  methods.  The  records  for  these  stations 
date  back  for  many  years  giving  the  Corps  an  excellent  understanding  of  how 
the  river  discharge  reacts  to  spring  thaws  and  large-scale  storm  events.  These 
records  provide  the  Corps  with  information  to  maintain  and  operate  the  navig¬ 
able  waterways  along  the  Mississippi.  As  technology  develops,  the  Corps  is 
always  eager  to  incorporate  new  technology  into  its  operating  plans,  but  not 
without  extensive  testing. 

The  Hydraulics  Laboratory  at  the  Waterways  Experiment  Station  was  con¬ 
tacted  by  the  Vicksburg  District  to  evaluate  the  present  discharge  measurement 
technique  against  a  new  technique  using  an  Acoustic  Doppler  Current  Profiler. 
A  testing  plan  was  developed  to  compare  these  two  techniques  over  two 
different  flow  conditions.  The  focus  of  this  report  is  to  present  the  data  from 
both  exercises,  describing  how  the  two  methods  calculate  discharge  through  a 
cross  section  and  identifying  possible  sources  of  error. 
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2  Acoustic  Doppler  Current 
Meter  Discharge 
Measurement  Technique 


An  Acoustic  Doppler  Current  Profiler  (ADCP)  is  an  instrument  that  uses 
acoustic  transducers  to  measure  the  movement  of  water  through  a  waterway. 
The  principles  of  operation  are  based  on  the  phenomena  of  doppler  shift  in 
signal  frequency  due  to  the  motion  of  a  second  frame  of  reference,  the  water. 
The  ADCP  sends  out  an  acoustic  signal  which  reflects  off  the  bottom.  Since 
the  water  and  the  ADCP  are  moving  independent  of  each  other,  the  ADCP 
measures  a  shift  in  the  signal  frequency  as  it  returns  through  the  water  column. 
Through  a  series  of  complicated  processing  algorithms,  the  frequency  shift  is 
related  to  the  individual  velocity  vector  components  of  the  water.  By  tracking 
the  boat’s  motion  the  ADCP  is  able  to  remove  the  boat’s  velocity  component 
from  the  water  velocity  measurement.  In  addition  to  a  velocity  measurement, 
the  ADCP  also  gives  a  bathymetry  profile  for  the  water  column  as  the  boat  is 
moving  along. 

There  are  several  different  types  of  ADCPs  which  differ  by  frequency  and 
size  of  transducers.  The  different  frequencies  allow  you  to  operate  in  different 
water  depth  ranges  with  varying  degrees  of  vertical  resolution.  Table  1  gives 
the  various  frequencies  available,  their  operating  ranges,  their  vertical 
resolution,  and  space  between  data  points  in  the  vertical.  The  different 
frequencies  allow  the  user  to  collect  data  over  a  wide  range  of  water  depths. 

An  ADCP  has  the  ability  to  incorporate  data  from  a  Global  Positioning 
System  (GPS)  in  the  data  collection  process  to  remove  the  boat’s  motion  from 
the  velocity  measurement.  The  user  has  an  option  of  referencing  the  bottom 
track  or  GPS  data  to  remove  the  boat’s  motion.  If  the  ADCP  cannot  track  the 
bottom,  it  cannot  calculate  a  discharge  for  the  cross  section.  The  bottom 
tracking  feature  provides  the  bathymetry  information  to  define  the  cross- 
sectional  area  for  the  discharge  calculation.  The  GPS  navigational  data  can 
provide  the  data  to  remove  the  boat’s  motion  from  the  velocity  measurement, 
but  an  external  data  set  from  a  fathometer  will  have  to  provide  the  data  for  the 
bottom  bathymetry.  The  two  data  sets  will  then  be  post-processed  to  make  a 
discharge  measurement.  This  will  be  the  only  option  available  for  discharge 
measurements  when  the  flows  are  too  high  for  bottom  tracking  to  be  effective. 
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The  physical  phenomena  that  causes  loss  of  bottom  tracking  is  movable  bed 
material.  When  the  water  velocities  increase  to  the  point  where  a  layer  of  sand 
or  mud  starts  to  move  across  the  bottom  in  sheets  or  waves,  then  the  ADCP 
has  trouble  distinguishing  between  the  channel  bottom  and  the  first  good  water 
layer  in  the  vertical.  There  are  several  ways  to  compensate  for  this  problem. 
One  method  is  to  increase  the  power  of  the  signal  to  punch  through  this  sheet 
or  wave  of  sediment  to  the  hard  bottom.  This  requires  using  a  lower- 
frequency  transducer  which  can  be  accomplished  by  switching  the  head  of  the 
transducer  or  using  a  completely  separate  ADCP.  Another  option  is  to  adjust 
some  of  the  operating  parameters  from  their  standard  default  values  to  improve 
bottom  tracking.  This  method  should  only  be  used  for  minor  bottom  tracking 
problems.  The  final  option  is  to  collect  a  second  data  set  of  the  bathymetry 
while  using  GPS  data  to  remove  the  boat’s  motion  from  the  velocity 
calculation.  In  a  post-processing  mode  the  two  data  sets  can  be  combined  to 
make  a  discharge  calculation. 

A  discharge  measurement  using  the  ADCP  is  accomplished  by  traversing 
the  river  or  stream  at  a  given  cross-section.  The  bathymetry  profile  coupled 
with  the  water  velocity  data  enables  the  ADCP  to  update  the  discharge 
calculation  as  the  boat  moves  across  the  waterway.  For  a  more  detailed 
explanation  of  how  the  TRANSECT  software  makes  the  discharge  calculation, 
Appendix  F  of  the  "User’s  Manual  for  the  RD  INSTRUMENTS  TRANSECT 
PROGRAM"*  has  been  included  as  Appendix  A  to  this  report. 


‘  RD  Instruments.  (1992).  "User’s  Manual  for  the  RD  INSTRUMENTS  TRANSECT 
PROGRAM,"  RD  Instruments,  San  Diego,  CA,  Appendix  F. 
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3  Current  Meter  Discharge 
Measurement  Technique 


The  Corp  presently  uses  the  current  meter  discharge  measurement 
technique^  on  established  discharge  ranges  on  large  rivers  and  streams.  It 
consists  of  summing  individual  discharge  calculations  across  the  discharge 
range.  The  partial  discharge  is  calculated  by  the  product  of  an  average 
velocity  measurement  and  an  average  area. 


Q  =  E(av) 


(Eq.  1) 


Q  represents  the  total  discharge  for  the  cross  section,  a  is  an  individual  partial 
cross-section  area  and  v  is  the  corresponding  average  velocity  normal  to  the 
partial  area. 

The  partial  area  is  that  of  a  rectangle,  with  a  width  one-half  way  between 
the  observation  point  before  and  after  the  current  observation  point  and  a  depth 
equal  to  that  of  the  current  observation  point  The  average  velocity  is  taken  at 
4/10-depth  of  the  water  column  for  a  predetermined  averaging  period. 

Figure  1  is  a  graphical  representation  of  the  river  cross  section.  The  cross 

section  is  defined  by  depths  at  locations  1,  2,  3,  4,  5 . n.  At  each  location 

the  velocities  are  measured  to  give  a  representative  average  for  the  entire  water 
column.  The  partial  discharge  for  any  partial  section  can  be  calculated  at  any 
location  x  by 


Q  =  V 

X 


(K  -  Vd)  .  (Vn  -  K) 


*  Thomas  J.  Buchanan  and  William  P.  Somers.  (1969).  "Discharge  measurements  at  Gaging 
Stations."  Techniques  of  water-resources  investigations  of  the  United  States  Geological  Survey, 

Book  3,  Chapter  A8,  U.S.  Government  Printing  Office,  Washington,  DC.  1-5. 
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Q  =  V 

X 


K  ♦  1)  ~  K  - 1) 


(Eq.2) 


where 


Q,  =  discharge  through  the  partial  section  x, 
y,  =  mean  velocity  at  location  jc, 

=  distance  from  initial  point  to  location, 
b(x.])  =  distance  from  the  initial  point  to  preceding  location, 
b(x*i)  =  distance  from  the  initial  point  to  next  location, 
dx  =  depth  of  water  at  location  x. 

The  discharge  through  partial  section  4  (heavily  outlined  in  Figure  1)  is 


b,  +  *3 
2 


When  X  is  at  the  beginning  or  end  of  the  cross  section,  then  equations  3  or  4 
are  used  to  calculate  the  partial  discharge. 


(Eq.  3) 


Q  =V 


b.  -  b. 


(» - 1) 


(Eq.  4) 


Figure  1  gives  an  example  of  a  possible  cross  section;  dj  is  zero  because 
the  depth  at  the  observation  point  1  is  zero.  When  the  cross-section  boundary 
is  a  vertical  line  at  the  edge  of  the  water  as  shown  at  location  n,  the  depth  is 
not  zero  and  the  velocity  at  the  end  section  may  or  may  not  be  equal  to  zero. 
Equations  3  and  4  are  used  whenever  there  is  water  only  on  one  side  of  an 
observation  point  such  as  piers,  abutments,  and  islands.  It  is  usually 
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necessary  to  estimate  the  velocity  at  an  end  section  because  it  normally  is 
impossible  to  measure  the  velocity  accurately  with  the  current  meter  close  to  a 
boundary.  The  summation  of  the  discharges  for  all  the  partial  sections  is  the 
total  discharge  of  the  stream. 


1,2,3, . « 

61,62,63, . 6n 

<[l,d2,d3, . dn 

Dashed  lines 


Observation  points 

Distance,  in  feet,  from  the  initial  point 
to  the  observation  point 

Depth  of  water,  in  feet,  at  the 
observation  point 

Boundary  of  the  partial  sections; 
one  heavily  outlined  discussed  in  text 


Figure  1 .  Graphical  representation  of  river  cross  section 
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4  Comparison  Plan  for  the 
Two  Discharge  Methods 


This  study  was  designed  to  compare  discharge  data  from  the  same 
discharge  ranges,  at  the  same  time  using  two  different  methods.  The  purpose 
of  this  effort  was  to  show  the  difference  and  precision  of  the  two  techniques. 
This  study  does  not  address  the  accuracy  of  the  two  methods  because  the  col¬ 
lection  sites  were  natural  sites  with  too  many  uncontrollable  factors. 

(Accuracy  refers  to  the  ability  to  measure  an  absolute  known  value  while  pre¬ 
cision  is  the  ability  to  repeat  the  measurement  regardless  of  the  accuracy.  In 
order  to  demonstrate  accuracy  of  the  two  methods,  a  controlled  flume  test  with 
a  known  discharge  would  have  to  be  conducted.)  The  difference  between  the 
two  methods  was  calculated  by  subtracting  the  mean  standard  method  from  the 
mean  ADCP  Method.  Then  the  percentage  that  the  difference  represents  for 
the  individual  means  was  calculated.  The  data  and  statistics  are  presented  in 
Tables  2-11. 

The  collection  effort  was  staged  during  two  different  flow  conditions, 
approximately  500,000  cfs  and  250,000  cfs.  We  chose  to  collect  data  for  these 
flows  to  see  if  there  were  any  significant  changes  in  the  difference  and  preci¬ 
sion  of  the  two  techniques.  The  boat  operating  conditions  were  different  for 
the  two  flow  conditions.  We  gained  insight  to  the  problems  associated  with 
making  discharge  measurements  on  the  river  under  various  flow  scenarios. 

By  comparing  the  two  methods  in  the  field,  we  have  been  able  to  identify 
possible  sources  of  error  in  both  techniques.  These  errors  need  to  be  further 
addressed  and  solutions  found  in  order  to  improve  both  techniques. 
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5  Data  from  ADCP  Discharge 
Measurement  Technique 


The  ADCP  data  were  collected  at  the  same  time  as  the  current  meter 
method  data.  The  files  were  processed  using  TRANSECT  software  in  the 
Playback  menu  option  in  order  to  take  advantage  of  the  estimating  capabilities 
of  the  software.  The  estimating  function  of  the  TRANSECT  software  can  only 
be  utilized  in  the  Playback  menu.  In  the  event  that  the  boat  is  unable  to  reach 
the  shoreline  because  of  shallow  water  or  some  other  obstruction,  the  user  can 
specify  the  distance  from  starting  and  ending  points  of  the  transect  line  to  the 
respective  shoreline  when  reprocessing.  The  software  uses  the  last  several 
good  ensembles  of  data  and  extrapolates  a  linear  fit  to  water’s  edge. 

Discharge  data  were  collected  at  four  discharge  ranges,  with  one  range 
being  repeated  on  the  fifth  day  of  collection.  This  sequence  of  field 
monitoring  was  repeated  for  a  mean  and  low-  flow  condition.  The  data  from 
the  two  collection  efforts  are  shown  in  Tables  2-11. 
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6  Data  from  the  Current 
Meter  Technique 


The  current  meter  data  were  collected  at  the  same  time  as  the  ADCP  data. 
The  contractor  under  Vicksburg  District  supervision  collected  three  data  sets  at 
each  range.  No  further  processing  was  done  to  these  data  after  the  contractor 
provided  the  boat  sheets  to  the  district.  The  discharge  data  collected  by  the 
contractor  are  presented  in  Tables  12-38.  The  total  values  for  discharge  were 
used  to  perform  simple  statistics  on  both  data  sets  which  are  presented  in 
Tables  2-11. 
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7  Sources  of  Error  for  Both 
Techniques 


Both  methods  have  proven  to  be  repeatable  in  discharge  measurements,  but 
what  causes  the  difference  between  the  two  methods?  There  are  two  types  of 
error  associated  with  each  method— user  error  and  sampling  error.  User  error  is 
the  only  type  of  error  that  can  be  addressed  with  the  present  data  set.  To 
identify  sampling  errors,  a  further  study  would  have  to  be  conducted  to 
quantify  method  accuracy. 

What  is  user  error?  User  error  is  any  type  of  error  associated  with  the 
setup  or  use  of  the  equipment. 

Standard  Method  User  Errors 

1.  Adjustment  of  the  bearings  on  the  Price  Meter 

2.  Boat  operation  while  collecting  data 

3.  Positioning  error  at  collection  points 

4.  Data  recording 

5.  Math  errors  in  making  the  discharge  calculation 

6.  Directional  errors 

ADCP  Method  User  Errors 

1.  Instrument  setup  before  operation 

2.  Boat  operation 

3.  Bottom  tracking 

Standard  Method 

The  Standard  Method  uses  a  Price  Meter  to  measure  water  flow  magnitude 
not  direction.  These  meters  require  maintenance  and  calibration  to  operate 
correctly.  The  bearing  on  the  rotating  cups  has  to  be  adjusted  and  kept  clean 
from  debris  in  order  to  turn  freely.  If  the  bearing  is  too  tight,  the  threshold 
velocity  will  be  higher  than  when  properly  adjusted;  and  if  the  bearing  is  too 
loose,  the  meter  will  wobble  and  turn  slower  than  it  should.  Bearing 
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adjustment  is  critical  for  proper  operation.  If  improper  adjustment  causes  the 
velocity  measurement  to  be  off  .1  foot  per  second  with  the  cross-sectional 
areas  measured  during  these  tests,  the  discharge  could  vary  +  10,000  cfs  for 
the  larger  cross  sections  and  +  5000  cfs  for  the  smaller  cross  sections. 

Boat  operation  is  another  significant  source  of  error  for  this  technique.  The 
boat  has  to  maintain  position  in  the  flow  while  the  timed  measurement  is  being 
made.  The  operator  uses  two  posts  on  each  side  of  the  river  to  line  up  with 
while  maintaining  position.  As  the  distance  from  the  shoreline  increases,  the 
error  in  this  type  of  system  becomes  larger.  The  operator  cannot  detect  small 
displacements  from  the  measurement  position  until  the  total  displacement 
exceeds  a  measurable  value.  With  all  of  these  factors  weighing  against  the 
operator,  it  is  almost  impossible  to  maintain  a  steady  platform  throughout  the 
measurement  period.  The  boat  is  either  drifting  slightly  ahead  or  behind  the 
point  of  data  collection.  This  drift  could  be  as  high  as  .1  foot  per  second  and 
go  undetected  by  the  operators  causing  a  difference  of  +  10,000  cfs  for  the 
larger  cross  sections  and  +  5,000  cfs  for  the  smaller  cross  sections. 

Positioning,  data  recording,  and  computational  errors  can  vary  based  on  the 
operators  attention  to  details.  For  very  careful  operators,  this  source  of  error 
could  be  minimized,  but  it  will  always  be  a  human  factor  which  is  variable. 

The  Standard  Method  has  no  way  to  account  for  velocity  directional 
variations  on  the  measured  cross  section.  The  bed  geometry  plays  a  significant 
part  in  establishing  flow  fields.  Large-  and  small-scale  eddies  move  along  the 
shoreline  because  of  the  bed  geometry.  It  was  noted  during  data  collection 
that  the  direction  of  flow  could  even  flow  upriver.  This  method  uses  the  entire 
flow  magnitude  to  calculate  discharge  rather  than  the  component  perpendicular 
to  the  cross  section.  On  some  sections  of  the  cross  section,  the  net  discharge 
for  that  small  section  should  be  subtracted  from  rather  than  added  to  the  total 
discharge.  On  large-scale  rivers,  like  the  Mississippi,  the  location  of  gaging 
ranges  is  very  critical  because  of  the  directional  limitations  of  the  Standard 
Method. 


ADCP  Method 

The  ADCP  Method  requires  an  experienced  operator  to  setup  and  use  the 
equipment.  Parameter  settings  can  greatly  affect  the  data  quality  from  the 
instrument.  Sensitivity  studies  on  how  much  each  parameter  affects  the 
instrument  operation  have  not  been  conducted.  With  a  proper  quality 
assurance  plan,  the  parameter  settings  can  be  defined  for  different  measurement 
scenarios. 

The  ADCP  Method  works  on  the  assumption  that  the  instrument  can  track 
the  boat’s  motion  over  the  bottom  while  recording  a  bathymetry  profile. 
TRANSECT  software  marks  the  bad  data  when  bottom  tracking  is  lost,  so  the 
operator  can  make  a  decision  about  rerunning  the  line.  One  section  of  data  in 
a  deep  swift  river  can  represent  a  significant  portion  of  the  discharge.  When 
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bottom  tracking  was  lost  during  the  field  tests,  the  total  discharge  could  be  10- 
20  percent  lower  than  the  mean  discharge  of  the  previous  measurements. 

Three  factors  affect  bottom  tracking—bed  geometry,  boat  speed,  and  bed¬ 
load  transport.  The  geometry  of  the  bed  can  cause  a  loss  in  bottom  tracking  if 
the  bottom  drops  off  or  rises  up  with  a  very  steep  angle.  The  only  recourse 
for  this  type  of  situation  is  to  slow  the  boat  down  on  the  approach  to  the  bed 
form  or  move  the  discharge  range  to  a  more  suitable  location.  Boat  speed 
affects  bottom  tracking,  but  there  is  not  an  optimum  speed  for  all  discharge 
measurements  sites.  The  one  rule  that  fits  the  best  for  this  problem  is  to  move 
across  the  cross  section  as  slow  as  possible  maintaining  navigation  control. 

The  data  are  recorded  on  a  set  time  interval  so  the  less  distance  traveled  in  that 
time  interval  means  more  data  points  in  a  cross  section.  More  data  points 
improve  the  total  discharge  measurement.  Bed-load  transport  is  the  hardest 
problem  to  overcome.  In  swift  rivers  where  sand  and  silt  travel  along  the 
bottom,  the  ADCP  has  trouble  identifying  the  bottom  signal.  The  operator  can 
adjust  some  of  the  operating  parameters  when  this  problem  occurs,  but  if  these 
adjustments  do  not  improve  tracking,  the  only  choices  are  to  use  a  lower- 
frequency  transducer  or  collect  bathymetry  data  with  a  separate  logging  system 
and  reference  navigational  GPS  data  to  account  for  the  boat’s  motion.  In  a 
post-processing  calculation,  the  two  data  sets  would  be  merged  to  calculate  a 
resulting  discharge. 
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Conclusions 


The  data  presented  in  Tables  2-11  show  that  the  discharge  measurements 
using  the  ADCP  Method  varied  less  about  the  mean  for  that  method.  In  the 
cases  where  only  one  or  two  discharge  measurements  were  made  using  the 
Standard  Method,  a  variance  for  the  Standard  Method  could  not  be  accurately 
calculated.  The  ADCP  Method  measurements  varied  a  maximum  of  1.25%  of 
the  mean  but  usually  the  variance  was  between  .5  and  .75  percent.  The 
Standard  Method  measurements  varied  a  maximum  of  2.0  percent  from  the 
mean  with  the  average  between  1.2  and  1.5  percent. 

A  significant  portion  of  the  difference  in  the  two  methods  can  be  attributed 
to  the  user  errors  stated  in  the  previous  section.  The  Mean  ADCP  Method 
appears  to  be  lower  than  the  Means  of  the  Standard  Method  in  most  field 
comparisons.  The  Means  of  the  ADCP  Method  were  higher  than  the  Means  of 
the  Standard  Method  when  the  boat  operator  drifted  with  the  flow  rather  than 
engaging  the  engines  during  the  collection  period.  This  was  noted  in  the  field 
notes  for  the  Chicot  and  Natchez  ranges.  On  both  days  that  data  were 
collected  on  these  ranges,  the  operator  of  the  boat  was  different  from  the 
operator  on  all  the  other  collection  days.  The  data  collected  using  the  ADCP 
Method  were  greater  than  the  Standard  Method  as  noted  in  Tables  5  and  8. 
Maintaining  a  steady  platform  for  the  Standard  Method  measurement  is  the 
cause  of  this  difference.  The  different  boat  operators  had  different  approaches 
to  maintaining  position  at  the  measurement  site.  One  operator  would  pull  up 
to  the  site  while  the  other  drifted  back  through  the  site.  The  drifting  boat 
removed  a  portion  of  the  velocity  measurement,  therefore  reducing  the  total 
discharge  calculation.  GPS  positioning  linked  with  a  navigation  package 
would  reduce  the  movement  of  the  boat  during  the  collection  period. 

The  solution  to  collecting  and  maintaining  an  accurate  data  set  for 
discharge  ranges  is  to  develop  Quality  Assurance  Plans  and  Standard  Operating 
Procedures.  These  plans  will  maintain  consistency  and  minimize  user  errors  in 
the  collection  process.  A  plan  of  action  also  needs  to  be  developed  for 
operation  in  high  flows  where  bottom  tracking  becomes  a  problem. 

Discharges  will  have  to  be  calculated  in  a  post-processing  mode  where  a 
bathymetry  data  set  is  coupled  with  the  ADCP  velocity  data.  The  motion  of 
the  boat  will  be  accounted  for  by  referencing  the  ADCP  data  to  GPS 
navigation  position  data. 
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Additional  software  that  needs  to  be  developed  will  address  the  problem  of 
loosing  one  ensemble  of  data  in  the  cross-section.  Sometimes  bed  geometry 
makes  it  impossible  to  collect  data  on  a  vertical  bank.  A  software  solution  to 
this  problem  would  use  the  data  on  either  side  of  the  lost  ensemble  to  fill  in 
the  lost  data.  This  would  address  bed  geometry  problems  and  allow  discharges 
to  be  corrected  in  a  post-processing  mode. 

The  ADCP  Method  should  enable  the  district  to  collect  more  consistent 
data  with  less  user  error  once  a  SOP  is  developed  for  the  discharge  program. 
The  additional  data  processing  options  for  the  increased  data  set  provided  with 
the  ADCP  Method  will  enable  the  district  to  make  better  engineering  decisions 
about  the  Mississippi  River  and  its  operation. 
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Table  1 

Types  of  ADCPs  with  Some  of  Their  Operational 
Characteristics 


Operating  Minimum  Bottom  Top 

Frequencies  of  Depth  Range  Bin  Size  Estimate  Estimate 

Transducers  Meters  Meters  Meters  Meters 


2400  KHz 

10 

0.12 

0.3 

0.3 

1200  KHz 

30 

0.25 

0.5 

0.5 

600  KHz 

100 

0.50 

1.0 

1.0 

300  KHz 

130 

1.00 

2.0 

2.0 

150  KHz 

230 

2.00 

4.0 

4.0 

75  KHz 


410 


4.00 


8.0 


8.0 


Table  2 

ADCP  Discharge  Data  on  the  Vicksburg  Discharge 


Range 


Transect 

Number 


Date 


Jul  11,  1994 


Jul  11,  1994 


Jul  11,  1994 


Jul  11,  1994 


Jul  11,  1994 


Jul  11,  1994 


Jul  11,  1994 


Jul  11,  1994  8 


Jul  11,  1994 


Jul  11,  1994  10 


Standard  deviation  of  measurements 


Mean  of  ADCP  method  measurements 


Length  of 
Transect  ft 


Course  Made 
Good  ft 


River  Discharge 
cu  ft/sec 


2687 

2625 

503922 

2688 

2633 

505702 

2686 

2639 

500915 

2754 

2644 

505119 

2805 

2667 

503656 

2757 

2660 

506428 

2776 

2671 

498509 

2754 

2690 

503268 

2741 

2658 

491939 

2725 

2657 

497501 

Percentage  of  the  mean  that  the  standard 
deviation  represents 


Discharge  Data  From  the  Standard  Collection  Methods 


19.49 


2654.4 


.734% 


3146.7 


502780.6 


.625% 


Cross-section 

Length 


Jul  11,  1994 

1 

Jul  11,  1994 

2 

Jul  11,  1994 

3 

River  Discharge 
cu  ft/sec 


Standard  deviation  of  measurements 


Mean  of  standard  method  measurements 


Percentage  of  the  mean  that  the  standard  deviation  represents 


Difference  of  mean  discharges  (ADCP  Method  -  Standard  Method) 


ADCP  Method 


Percentage  of  the  method  means  that  the 
difference  represents 


*  Some  data  were  lost  on  the  cross  section 


Standard 


526870 


535245 


522146 


6633.76 


528087 


1 .26% 


-25307 


5.03% 


4.79% 


Table  3 

ADCP  Discharge  Data  on  the  Natchez  Discharge  Range 

Date 

Transect 

Number 

Length  of 
Transect  ft 

Course  Made 
Good  ft 

River  Discharge 
cu  ft/sec 

Jul  12,  1994 

1 

2842 

2794 

512346 

Jul  12,  1994 

2 

2884 

2807 

508610 

Jul  12,  1994 

3 

2851 

2798 

507101 

Jul  12,  1994 

4 

2849 

2790 

517517 

Jul  12,  1994 

5 

2843 

2800 

518188 

Jul  12,  1994 

6 

2848 

2789 

510710 

Jul  12,  1994 

7 

2861 

2784 

501887 

Jul  12,  1994 

8 

2876 

2798 

507623 

Jul  12,  1994 

9 

2931 

2806 

504965 

Jul  12,  1994 

10 

2847 

2786 

511322 

Standard  deviation  of  measurements 

7.96 

5147.3 

Mean  of  ADCP  method  measurements 

2795,2 

510026.9 

Percentage  of  the  mean  that  the  standard 
deviation  represents 

0.27% 

1 .00% 

Discharge  Data  From  the  Standard  Collection  Methods 

Cross-section 

Length 

River  Discharge 
cu  ft/sec 

Jul  12  1994 

1 

2874 

534245 

Jul  12,  1994 

2 

2874 

515197 

Jul  12,  1994 

3 

2874 

534245 

Standard  deviation  of  measurements 

10997 

Mean  of  standard  method  measurements 

527895.7 

Percentage  of  the  mean  that  the  standard  deviation  represents 

2.08% 

Difference  of  mean  discharges  (ADCP  Method  -  Standard  Method) 

-17868.8 

Percentage  of  the  method  means  that  the 
difference  represents 

ADCP  Method 

3.50% 

Standard 


3.38% 


Table  4 

ADCP  Discharge  Data  on  the  Grand  Gulf-Kempe  Bend  Discharge 
Range 


Date 

Transect 

Number 

Length  of 
Transect  ft 

Course  Made 
Good  ft 

River  Discharge 
cu  ft/sec 

Jul  13,  1994 

1 

5367 

5279 

464691 

Jul  13,  1994 

2 

5457 

5299 

467459 

Jul  13,  1994 

3 

5411 

5311 

462104 

Jul  13,  1994 

4 

5408 

5303 

459282 

Jul  13,  1994 

5 

5511 

5316 

458964 

Jul  13,  1994 

6 

5441 

5321 

464559 

Jul  13,  1994 

7 

5482 

5301 

451080 

Jul  13,  1994 

8 

5464 

5304 

461161 

Jul  13,  1994 

9 

5380 

5285 

456593 

Jul  13,  1994 

10 

5499 

5304 

443645 

Standard  deviation  of  measurements 

31.15 

3585.0 

Mean  of  ADCP  method  measurements 

5292.3 

461851.6 

Percentage  of  the  mean  that  the  standard 
deviation  represents 

0.67% 

0.78% 

Discharge  Data  From  the  Standard  Collection  Methods 


Cross-section 

Length 

River  Discharge 
cu  ft/sec 

5388 

494143 

5388 

504693 

5388 

500907 

Jul  13,  1994 

1 

Jul  13,  1994 

2 

Jul  13,  1994 

3 

Standard  deviation  of  measurements 


Mean  of  standard  method  measurements 


Percentage  of  the  mean  that  the  standard  deviation  represents 


Difference  of  mean  discharges  (ADCP  Method  -  Standard  Method) 


ADCP  Method 


Percentage  of  the  method  means  that  the 
difference  represents 


Standard 


5344.2 


499914.3 


1.07% 


-38062.7 


8.24% 


7.61% 


Some  data  were  lost  on  the  cross  section 


Table  5 

ADCP  Discharge  Data  on  the  Chicot  Discharge  Range 


Date 

Transect 

Number 

Length  of 
Transect  ft 

Course  Made 
Good  ft 

River  Discharge 
cu  ft/sec 

Jul  14,  1994 

1 

2945 

2911 

429212 

Jul  14,  1994 

2 

2957 

2903 

425371 

Jul  14,  1994 

3 

2922 

2822 

426260 

Jul  14,  1994 

4 

2908 

2847 

428157 

Jul  14,  1994 

5 

2896 

2838 

426807 

Jul  14,  1994 

6 

2931 

2858 

428859 

Jul  14,  1994 

7 

2933 

2861 

432886 

Jul  14,  1994 

8 

2927 

2851 

439166 

Jul  14,  1994 

9 

2903 

2850 

426793 

Jul  14,  1994 

10 

2906 

2856 

431149 

Standard  deviation  of  measurements 

27.38 

4110.0 

Mean  of  ADCP  method  measurements 

2859.7 

429466 

Percentage  of  the  mean  that  the  standard 
deviation  represents 

0.96% 

0.96% 

Discharge  Data  From  the  Standard  Collection  Methods 


Cross-section 

Length 

River  Discharge 
cu  ft/sec 

2914 

409097 

2914 

415038 

2914 

410178 

Jul  14,  1994 

1 

Jul  14,  1994 

2 

Jul  14,  1994 

3 

Standard  deviation  of  measurements  3164.48 


Mean  of  standard  method  measurements  41 1437.7 


Percentage  of  the  mean  that  the  standard  deviation  represents  0.77% 


Difference  of  mean  discharges  (ADCP  Method  -  Standard  Method)  18028.3 


Percentage  of  the  method  means  that  the  ADCP  Method  4.20% 

difference  represents 


Standard 


4.38% 


Table  6 

ADCP  Discharge  Data  on  the  Vicksburg  Discharge  Range 


Date 

Transect 

Number 

Length  of 
Transect  ft 

Course  Made 
Good  ft 

River  Discharge 
cu  ft/sec 

Jul  15,  1994 

1 

2581 

2492 

456883 

Jul  15,  1994 

2 

2577 

2471 

444848 

Jul  15,  1994 

3 

2622 

2513 

449096 

Jul  15,  1994 

4 

2592 

2504 

455911 

Jul  15,  1994 

5 

2612 

2524 

455221 

Jul  15,  1994 

6 

2616 

2528 

452352 

Jul  15,  1994 

7 

2607 

2513 

453573 

Jul  15,  1994 

8 

2578 

2507 

448193 

Jul  15,  1994 

9 

2603 

2505 

451752 

Jul  15,  1994 

10 

2600 

2513 

451988 

Standard  deviation  of  measurements 

16.23 

3745.3 

Mean  of  ADCP  method  measurements 

2507 

451981.7 

Percentage  of  the  mean  that  the  standard 
deviation  represents 

0.65% 

0.82% 

Discharge  Data  From  the  Standard  Collection  Methods 


Jul  15,  1994 

1 

Jul  15,  1994 

2 

Jul  15,  1994 

3 

Cross-section 

River  Discharge 

Length 

cu  ft/sec 

Standard  deviation  of  measurements 


Mean  of  standard  method  measurements 


Percentage  of  the  mean  that  the  standard  deviation  represents 


Difference  of  mean  discharges  (ADCP  Method  -  Standard  Method) 


ADCP  Method 


Percentage  of  the  method  means  that  the 
difference  represents 


Standard 


471127 


469147 


453150 


2001.5 


471141.3 


0.42% 


-19159.6 


4.23% 


4.06% 


Table  7 

ADCP  Discharge  Data  on  the  Vicksburg  Discharge  Range 

Date 

Transect 

Number 

Length  of 
Transect  ft 

Course  Made 
Good  ft 

River  Discharge 
cu  ft/sec 

Sep  19,  1994 

1 

2048 

1998 

254494 

Sep  19, 1994 

2 

2068 

2025 

260393 

Sep  19, 1994 

3 

2117 

2097 

256447 

Sep  19, 1994 

4 

2101 

2072 

256003 

Sep  19, 1994 

5 

2075 

2043 

259188 

Sep  19,  1994 

6 

2057 

2036 

257120 

Sep  19,  1994 

7 

2115 

2041 

261795 

Sep  19,  1994 

8 

2147 

2093 

258083 

Sep  19,  1994 

9 

2119 

2081 

257260 

Sep  19,  1994 

10 

2104 

2067 

256611 

Standard  deviation  of  measurements 

31.91 

2180.9 

Mean  of  ADCP  method  measurements 

2055.3 

257739.4 

Percentage  of  the  mean  that  the  standard 
deviation  represents 

1.55% 

0.84% 

Standard  deviation  of  measurements 

1180.86 

Mean  of  standard  method  measurements 

282435 

Percentage  of  the  mean  that  the  standard  deviation  represents 

0.41  % 

Difference  of  mean  discharges  (ADCP  Method  -  Standard  Method) 

-24695.6 

Percentage  of  the  method  means  that  the 
difference  represents 

ADCP  Method 

9.58% 

Standard 


8.74% 


Table  8 

ADCP  Discharge  Data  on  the  Natchez  Discharge  Range 


Date 

Transect 

Number 

Length  of 
Transect  ft 

Course  Made 
Good  ft 

River  Discharge 
cu  ft/sec 

Sep  20,  1994 

1 

2760 

2723 

272144 

Sep  20,  1994 

2 

2771 

2725 

271658 

Sep  20,  1994 

3 

2743 

2705 

270884 

Sep  20,  1994 

4 

2749 

2729 

272353 

Sep  20,  1994 

5 

2737 

2703 

270579 

Sep  20, 1994 

6 

2765 

2724 

271808 

Sep  20,  1994 

7 

2759 

2718 

272126 

Sep  20,  1994 

8 

2734 

2704 

270611 

Sep  20,  1994 

9 

2756 

2727 

272785 

Sep  20,  1994 

10 

2773 

2727 

271258 

Standard  deviation  of  measurements  10.43 


Mean  of  ADCP  method  measurements  2718.5 


Percentage  of  the  mean  that  the  standard  0.38% 

deviation  represents 


Discharge  Data  From  the  Standard  Collection  Methods 


762.77 


271620.6 


0.28% 


Sep  20,  1994 

1 

Sep  20,  1 994 

2 

Sep  20,  1994 

3 

Cross-section 

River  Discharge 

Length 

cu  ft/sec 

2788 

270007 

2788 

267787 

Standard  deviation  of  measurements 


Mean  of  standard  method  measurements 


Percentage  of  the  mean  that  the  standard  deviation  represents 


Difference  of  mean  discharges  (ADCP  Method  -  Standard  Method) 


Percentage  of  the  method  means  that  the 
difference  represents 


1569.77 


268897 


0.58% 


2723.6 


ADCP  Method 

1.00 

Standard 

1 .01  % 

Table  9 

ADCP  Discharge  Data  on  the  Grand  Gulf-Kempe  Bend  Discharge 
Range 


Date 

Transect 

Number 

Length  of 
Transect  ft 

Course  Made 
Good  ft 

River  Discharge 
cu  ft/sec 

Sep  21,  1994 

1 

2083 

2041 

240468 

Sep  21,  1994 

2 

2138 

2073 

246464 

Sep  21,  1994 

3 

2104 

2055 

240715 

Sep  21,  1994 

4 

2138 

2069 

238617 

Sep  21,  1994 

5 

2119 

2078 

243221 

Sep  21,  1994 

6 

2148 

2096 

246042 

Sep  21,  1994 

7 

2084 

2044 

239550 

Sep  21,  1994 

8 

2128 

2076 

245615 

Sep  21,  1994 

9 

2115 

2085 

241009 

Sep  21,  1994 

10 

1932 

1889 

243203 

Standard  deviation  of  measurements 

10.43 

2834.2 

Mean  of  ADCP  method  measurements 

2718.5 

242490.4 

Percentage  of  the  mean  that  the  standard 
deviation  represents 

0.38% 

1.17% 

Discharge  Data  From  the  Standard  Collection  Methods 


Sep  21,  1994 

1 

Sep  21,  1994 

2 

Sep  21,  1994 

3 

Cross-section 

Length 


River  Discharge 
cu  ft/sec 


2220 

262307 

2220 

264399 

Standard  deviation  of  measurements 


Mean  of  standard  method  measurements 


Percentage  of  the  mean  that  the  standard  deviation  represents 


Difference  of  mean  discharges  (ADCP  Method  -  Standard  Method) 


ADCP  Method 


Percentage  of  the  method  means  that  the 
difference  represents 


Standard 


1479.26 


263353 


0.56% 


-20862.6 


8.60% 


7.92% 


Table  10 

ADCP  Discharge  Data  on  the  Chicot  Discharge  Range 


Transect 

Number 


Length  of 
Transect  ft 


Date 


Sep  29,  1994  1 


Sep  29, 1994 


Sep  29,  1994 


Sep  29,  1994 


Sep  29,  1994 


Sep  29,  1994 


Sep  29,  1994 


Sep  29,  1994  8 


Sep  29, 1994 


Sep  29,  1994  10 


Standard  deviation  of  measurements 


Mean  of  ADCP  method  measurements 


Percentage  of  the  mean  that  the  standard 
deviation  represents  _ 


Discharge  Data  From  the  Standard  Collection  Methods 


Course  Made 
Good  ft 

River  Discharge 
cu  ft/sec 

2794 

236675 

2783 

236790 

2773 

240504 

2791 

236753 

2783 

236314 

2791 

235972 

2789 

237100 

2667 

236470 

2776 

236203 

2768 

237359 

37.7 

1293.96 

2771.5 

237014 

1 .36% 

0.54% 

Cross-section 

Length 


River  Discharge 
cu  ft/sec 


Sep  29,  1994  1 


Sep  29,  1994  2 


Sep  29,  1994  3 


Standard  deviation  of  measurements 


Mean  of  standard  method  measurements _ _ 


Percentage  of  the  mean  that  the  standard  deviation  represents _ 


Difference  of  mean  discharges  (ADCP  Method  -  Standard  Method) 


Percentage  of  the  method  means  that  the  ADCP  Method 

difference  represents 


248438 


243870 


246531 


2294.34 


246279.7 


0.93% 


-9265.7 


3.91% 


Standard 


3.76% 


Table  1 1 

ADCP  Discharge  Data  on  the  Vicksburg  Discharge  Range 


Transect 

Number 


Length  of 
Transect  ft 


Course  Made 
Good  ft 


Sep 

23,  1994 

Sep 

23,  1994 

Sep 

23,  1994 

Sep 

23, 1994 

Sep 

23,  1994 

Sep 

23,  1994 

Sep 

23, 1994 

Sep 

23,  1994 

Sep 

23,  1994 

Sep 

23,  1994 

River  Discharge 
cu  ft/sec 


2088 

2051 

252845 

2094 

2035 

252234 

2088 

2055 

245334 

2076 

2054 

240640 

2103 

2041 

245002 

2073 

2036 

248952 

2081 

2050 

245581 

2105 

2058 

252723 

2090 

2053 

247079 

2103 

2067 

246096 

Standard  deviation  of  measurements  10.03 


Mean  of  ADCP  method  measurements  2050 


Percentage  of  the  mean  that  the  standard  0.49% 

deviation  represents 


Discharge  Data  From  the  Standard  Collection  Methods 


3342.7 


248427.3 


1.34% 


Cross-section 

Length 

River  Discharge 
cu  ft/sec 

2202 

253450 

2202 

254595 

2202 

254994 

Sep  29,  1994 

1 

Sep  29,  1994 

2 

Sep  29,  1994 

3 

Standard  deviation  of  measurements 


Mean  of  standard  method  measurements 


Percentage  of  the  mean  that  the  standard  deviation  represents 


Difference  of  mean  discharges  (ADCP  Method  -  Standard  Method) 


ADCP  Method 


Percentage  of  the  method  means  that  the 
difference  represents 


Standard 


801 .47 


254346.3 


0.32% 


-5919 


2.38% 


2.32% 
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Table  13 

Range  Gage  435.41 ,  Vicksburg  Discharges 


Notes;  All  distances  from  AB  at  right  bank: 

AB  to  RWE  =  434  ft 
AC  to  LWE  =  277  ft 

Stream  discharge  measurements,  1 1  July  1994,  mean  time  1020,  width  2,710  ft,  area  1 18,000  ft^  mean  velocity  4.53  fps,  maximum  velocity  6.57  fps, 
minimum  velocity  2.29  fps,  discharge  535,000  ft^/min,  gage  height  20.93,  A  gage  -0.08,  method  coefficient  0.4,  measurement  interval  300  sec. 


Table  15 

Range  Gage  362.34,  Natchez  Discharges 


Notes;  All  distances  from  AB  at  right  bank; 

ABto  RWE  =  104  ft 
AC  to  LWE  =  71  ft 

Stream  discharge  measurements,  12  July  1994,  mean  time  1045,  width  2,870  ft,  area  143,000  ft^  mean  veiocity  3.60  fps,  maximum  velocity  4.93  fps, 
minimum  velocity  1 .81  fps,  discharge  515,000  ft^/min,  gage  height  27.63,  A  gage  -0.08,  method  coefficient  0.4,  measurement  interval  300  sec. 


Table  21 

Range  Gage  565.9,  Chicot  Discharges 


Table  22  (Concluded) 


00 

o> 


o 

o> 


Table  23 

Range  Gage  565.9,  Chicot  Discharges 


Notes;  All  distances  from  AB  at  right  bank: 

Stream  discharge  measurements,  14  July  1994,  mean  time _ ,  width  2,910  ft,  area  130,000  ft^  mean  velocity  3.19  fps,  maximum  velocity  4.51  fps, 

minimum  velocity  0.40  fps,  discharge  415,000  ft^/min,  gage  height  13.01,  A  gage  +0.06,  method  coefficient  0.4,  measurement  interval  300  sec. 


Table  23  (Concluded) 


Table  34 

Range  Gage  565.9,  Chicot  Discharges 


Notes:  All  distances  from  AB  at  right  bank: 

Stream  discharge  measurements,  29  Sept  1994,  mean  time  1030,  width  2,830  ft,  area  109,000  ft^  mean  velocity  2.24  fps,  maximum  velocity  3.31  fps, 
minimum  velocity  0.56  fps,  discharge  244,000  ft^/min,  gage  height  2.98,  A  gage  -0.04,  method  coefficient  0.4,  measurement  interval  300  sec. 


Table  35 

Range  Gage  565.9,  Chicot  Discharges 


Notes;  All  distances  from  AB  at  right  bank: 

Stream  discharge  measurements,  29  Sept  1994,  mean  time  1030,  width  2,830  ft,  area  109,000  fr  mean  velocity  2.27  fps,  maximum  velocity  3.13  fps, 
minimum  velocity  0.53  fps,  discharge  247,000  ft^/min,  gage  height  2.98,  A  gage  -0.04,  method  coefficient  0.4,  measurement  intervai  300  sec. _ 


Table  37 

Range  Gage  435.41 ,  Vicksburg  Discharges 


Notes:  All  distances  from  AB  at  right  bank: 

AB  to  RWE  =  870  ft 
AC  to  LWE  =  349  ft 

Stream  discharge  measurements,  23  Sept  1994,  mean  time  1020,  width  2,200  ft,  area  88,500  ft^  mean  velocity  2.88  fps,  maximum  velocity  4  50  fos 
minimum  velocity  0.72  fps,  discharge  255,000  ff^/min,  gage  height  8.33,  A  gage  -0.02,  method  coefficient  0.4,  measurement  interval  300  sec. 


Appendix  A 

ADCP  Discharge  Calcuiation 
Methods 


This  appendix  is  a  technical  description  of  how  the  ADCP  collects  and 
calculates  discharge  data.  This  appendix  has  been  published  as  Appendix  F  of 
"User’s  Manual  for  the  RD  INSTRUMENTS  TRANSECT  PROGRAM"'  and 
is  reprinted  with  permission  of  RD  Instruments.  The  appendix  includes  two 
change  releases.  On  a  change  page,  a  vertical  bar  in  the  left  margin  shows  the 
part  of  text  affected  by  the  change.  The  last  line  (footer)  of  each  page  lists  the 
effective  date  of  the  change. 


‘  RD  Instruments.  (1992).  "User’s  Manual  for  the  RD  INSTRUMENTS  TRANSECT 
PROGRAM,"  RD  Instruments,  San  Diego,  CA,  Appendix  F. 
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APPENDIX  F 

FORMULAS  AND  CALCULATIONS 


F-I.  DISCHARGE  CALCULATIONS 

This  seccion  explains  how  TRANSECT  calculates  discharge.  For  the  data  mea¬ 
sured  by  the  ADCP,  TRANSECT  calculates  discharge  (middle  water  layer  -  MidQ) . 
For  the  unmeasured  parts  of  the  profile  (top  water  layer  -  TopQ,  bottom  water 
layer  -  BtmQ)  TRANSECT  estimates  the  discharge.  TRANSECT  accumulates  these 
values  over  the  entire  transect  (or  subsection  of  a  transect  if  SUBSECTIONING  IS 
NOW  •  ON  in  the  PLAYBACK  menu) .  When  <alt-u>  unmeasured  is  active ,  TRANSECT  displays 
the  accumulated  values  on  the  right  side  of  the  screen.  The  total  discharge 
(2Q)  is  the  summation  of  discharge  in  the  top,  middle,  and  bottom  layers. 
TRANSECT  also  can  estimate  near-shore  discharge  and  included  in  the  total  dis¬ 
charge  in  the  playback  menu. 

Discharge  is  the  total  volume  of  water  flowing  through  a  cross-section  of 
water  per  unit  of  time.  TRANSECT  computes  this  total  volume  discharge  (SQ) 
for  each  ADCP  ensemble  and  displays  the  estimate  in  the  ACQUIRE  and  PLAYBACK 
modes.  An  ADCP  measures  profiles  of  water-current  velocity  relative  to  the 
vessel.  The  ADCP  also  measures  the  velocity  of  the  vessel  relative  to  the 
bottom  and  depth  to  the  bottom  for  each  ADCP  beam.  Computation  of  discharge 
depends  only  on  these  data.  We  do  not  need  to  know  compass  heading  or  vessel 
location.  Furthermore,  the  transect  can  be  an  arbitrary  curve  (see  Figure 
F-1)  as  long  as  it  starts  near  one  side  of  a  channel  and  ends  near  the  other. 

The  uncertainty  in  the  discharge  estimate  arises  from  random  errors,  biases, 
and  missed  data  (near  the  surface  and  bottom,  and  near  the  sides  of  a  chan¬ 
nel)  .  TRANSECT  can  extrapolate  near-shore  discharge  (near  the  channel  sides) 
during  playback.  The  algorithm  for  estimating  discharge  is  adopted  from 
Simpson  and  Oltmann  (1990),  and  Gordon  (1989). 

There  are  two  methods  available  in  TRANSECT  for  estimating  the  discharge  in 
the  unmeasured  parts  of  the  profile.  You  can  use  either  the  CONSTANT  method 
(i.e.  ,  "straight  up  and  down")  or  a  POWER  law  method.  If  you  select  a  POWER 
method  for  either  the  top  or  bottom  unmeasured  part  of  the  profile,  you  can 
select  the  exponent  of  the  power  law  based  on  flow  conditions  and  the  rough¬ 
ness  of  the  channel  bed.  For  more  information  on  the  theory  of  power  laws  for 
flow  resistance  see  Chen  (1991). 

For  this  version  of  TRANSECT,  bottom-track  velocity  data  from  the  ADCP  must  be 
valid  for  the  moving-vessel  discharge  calculation  to  be  correct.  The  dis¬ 
charge  profile  plot  in  the  ACQUIRE  and  PLAYBACK  modes  displays  the  depth  to 
which  TRANSECT  is  measuring  discharge  for  the  current  ADCP  ensemble.  The 
TABULAR  submenu  shows  the  discharge  value  and  its  corresponding  depth. 

We  begin  with  a  discussion  of  the  moving-vessel  method  of  determining  dis¬ 
charge  using  an  ADCP.  We  then  discuss  how  TRANSECT  implements  this  method  and 
determines  the  size  of  the  top,  bottom,  and  middle  water  layers.  Table  F-1 
lists  several  terms  used  throughout  our  discussion.  The  explanations  present¬ 
ed  here  assume  the  ADCP  is  looking  down,  so  first  implies  shallowest  and  last 
implies  deepest. 
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Table  F-1.  Discharge  Calculation  Terms 


Term 

Definition 

at  1 

Time  difference  between  successive  ensembles. 

BCmQ 

Estimated  discharge  for  the  unmeasured  data  at  the  bottom  of  the  profile. 

*D. 

Depth  cell  (bin)  size,  or  the  length  of  the  range  gate.  Each  range 
gate  corresponds  to  a  depth  cell.  The  ADCP  constructs  a  profile  from  a 
series  of  range  gates  or  depth  cells. 

^ADCP 

Depth  of  the  ADCP  transducer  face  from  the  water  surface. 

*Di, 

Blank  beyond  transmit. 

Dlc 

Depth  of  the  last  good  bin  (i.e.,  the  last  bin  in  middle  water  layer). 

^LCmax 

Depth  of  the  last  possible  good  bin. 

Depth  to  the  bottom  (i.e.,  the  channel  bed). 

D.V. 

Averaged  measured  depth  from  the  ADCP  beams  (not  including  D^^scp)  • 

*D„ 

Lag  between  transmit  pulses  or  correlation  lag. 

*Dp 

Transmit  pulse  length.  Dp  is  the  length  of  a  single  transmit  pulse. 

Broadband  systems  transmit  two  or  more  pulses.  If  the  pulse  is  coded, 
this  is  the  total  length  of  one  coded  pulse,  not  the  length  of  a  single 
element. 

*Dtop 

Depth  of  the  center  of  the  first  bin. 

MidQ 

Discharge  computed  for  the  middle  of  the  profile  from  ADCP  velocity  and 
bottom- track  data. 

*Np 

Number  of  transmit  pulses . 

TopQ 

Estimated  discharge  for  the  unmeasured  data  at  the  top  of  the  profile. 

X 

Cross-product  computed  from  ADCP  data  (see  i[F-l-2) 

X' 

Synthetic  cross-product  computed  from  the  power-curve  fit 

•NOTE:  THESE  DATA  ARE  OIRECTLT  OUTPUT  BY  THE  ADCP. 

F-2 

TRANSECT  USER’S  MANUAL  (B0)  -  OCTOBER  1992  (CHANGE  2) 

Appendix  A  ADCP  Discharge  Calculation  Methods 


A3 


APPENDIX  F  -  FORMULAS  AND  CALCUUTIONS 


F-1.1.  DETERMINING  MOVING-VESSEL  DISCHARGE  AND  THE  CROSS-PRODUCT.  The 

moving-vessel  method  for  measuring  total  discharge  using  an  ADCP  is  computed 
by  transecting  a  channel  from  bank  to  bank  and  accumulating  the  discharge  for 
each  ensemble  (see  Figure  F-1).  The  general  equation  for  determining  total 
channel  discharge  (Qt)  through  an  arbitrary  surface  (s)  is  (Simpson  and 
Oltmann,  1990;  Gordon,  1989):. 


II  Vf  •  n 


(Eq.l) 


where:  ds  -  Differential  area. 

Vf  -  Mean  water  velocity  vector. 

n  -  Unit  vector  normal  to  transect  path  at  a  general  point. 

For  a  vessel  moving  across  a  channel,  the  area  of  s  is  determined  by  the  ver¬ 
tical  surface  beneath  the  transect  path.  Because  the  ADCP  measures  both  ves¬ 
sel  motion  (ADCP  bottom- track  velocity)  and  water  velocity,  we  can  re-order 
the  integral  above  (Eq.l)  in  the  following  form. 

Let:  dz  -  Differential  depth 

dt  -  Differential  time 
Vj,  -  Mean  vessel  velocity  vector 
T  -  Total  transect  time 
d  -  Total  depth  (Dt„„i) 

k  -  Unit  vector  in  the  vertical  direction 


ds  -  Vk  dz  dt 


Vf.n  -  |Vf|  sin(«l) 


(Eq.2) 

(Eq.3) 


XU 

II  V£-n  ds  -  II  IVfl  IV, 


bl  sin(^)  dz  dt 


(Eq.4) 


(Vf  X  Vb)  •  k  dz  dt 


(Eq.5) 


Converting  the  cross-product  x  V^) 
nates,  we  get: 


k  into  rectangular  (vessel)  coordi- 


where : 


(Vf  XV,)  .  k  -  F.B^  -  F,B,.  (Eq.6) 

where:  F„  -  Cross -component  of  the  mean  water  velocity  vector. 

Fy  -  Fore/aft  component  of  the  mean  water  velocity  vector. 

B„  -  Cross -component  of  the  mean  vessel  velocity  vector. 

By  -  Fore/aft  component  of  the  mean  vessel  velocity  vector. 

The  values  for  F,,  Fy,  B^,  By  are  rotated  to  earth  coordinates  before  computing 
the  cross-product  above  (which  does  not  affect  the  value  of  the  cross- 
product).  Note  that  TRANSECT  converts  these  values  to  earth  coordinates  be¬ 
fore  any  processing  is  done.  Throughout  the  remainder  of  this  appendix,  we 
will  assign  the  symbol  X  to  represent  the  cross-product. 
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F-1.2.  METHODS  FOR  ESTIMATING  DISCHARGE  IN  THE  UNMEASURED  TOP/BOHOM  PARTS  OF 
THE  VELOCITY  PROFILE.  There  are  two  methods  available  in  TRANSECT  to  estimate 
discharge  in  the  unmeasured  top/bottom  parts  of  the  velocity  profile  based  on 
the  TOP/BOTTOH  DISCHARGE  ESTIMATE  settings  in  the  CFG  file  (set  in  the  CALIBRATION 
menu) .  The  two  methods  are  CONSTANT  and  POWER. 

Constant  Method.  This  is  the  simplest  method  of  estimating  the  discharge  in 
the  unmeasured  parts  of  the  profile.  However,  this  method  does  not  follow 
accepted  hydraulic  descriptions  of  the  vertical  distribution  of  horizontal 
water  velocities  in  open  channels,  particularly  in  the  bottom  water  layer  near 
the  bottom  boundary  where  the  velocity  decreases  to  zero.  Simpson  and  Oltmann 
(1990),  and  Gordon  (1989),  discuss  this  method.  This  "straight  up  and  down" 
method  extrapolates  the  cross-product  from  the  first  good  bin  in  the  profile 
to  the  top  (TopQ)  or  from  the  last  good  bin  to  the  bottom  (BtmQ) .  See  ^F-1.6 
and  ^F-1.7  for  details  on  how  TRANSECT  does  this. 


Power  Method.  Chen  (1991)  discusses  the  theory  of  power  laws  for  flow  resis 
tance.  Simpson  and  Oltmann  (1990)  discuss  Chen's  power  law  equivalent  of 
Manning's  formula  for  open  channels  (with  b  -  1/6). 


u/u’  -  9.5  (Z/Zo)*" 

(Eq.7) 

where:  z  -  Distance  to  the  channel  bed. 

u  -  Velocity  at  distance  z  from  bed. 
u*  -  Shear  velocity. 

Zo  -  Bottom  roughness  height, 
b  -  Exponent  (1/6). 

Combining  terms  not  from  the  ADCP  we  get: 

u  -  (9.5u7zo’*)z'’- 

(Eq.8) 

If  we  let  a'  -  (9.5u*/Zo’),  then: 

u  -  a'  z*" . 

(Eq.9) 

For  the  moving  boat  discharge  calculation,  the  cross' 
from  the  ADCP  velocities  (Eq.  6),  replacing  u  in  the 

■product  (X)  is  computed 
power  law  of  Eq.  9: 

X  -  a"z‘’. 

(Eq.lO) 

For  each  depth  cell,  the  distance  from  the  channel  bed  z  and  the  cross-product 
X  are  computed  by  TRANSECT.  The  next  step  is  to  solve  for  the  unknown  a"  in 
the  power  law  of  Eq.  10.  TRANSECT  solves  for  a"  by  setting  the  antiderivative 
of  the  power  law  equal  to  the  cross-product  values  integrated  over  the  valid 
depth  range  of  the  profile. 

Let:  Zj  -  Distance  from  the  channel  bed  to  last  (deepest)  good  bin 

Z2  -  Distance  from  the  channel  bed  to  first  (shallowest)  good  bin 

The  values  for  Zi  and  Z2  define  the  region  of  the  profile  with  valid  ADCP 
data  (see  Figure  F-2) .  This  is  referred  to  as  the  middle  layer  of  the  pro¬ 
file.  Section  F-1.4  discusses  how  TRANSECT  determines  Zj  and  Z2. 
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First,  we  integrate  the  ADCP  data  over  the  valid  depth  range. 

Let:  "  1 

Zi 

where:  dz  -  depth  cell  size  (D,) . 


(Eq.ll) 


therefore ; 


li  -  D,  E  X, 


(Eq.l2) 


Second,  use  the  antiderivative  of  the  power  law  of  Eq.  10  to  integrate  the 
power  law  over  the  middle  water  layer: 


f(z)  -  a^'z*”  =»  F(z)  - 


(Eq.l3) 


Rz  -  f(z)  dz 


(Eq.l4) 


(Eq.15) 


a" 


(Eq.l6) 


Equating  the  integrals  we  solve  for  a"; 


Ri  “  Rj 


(Eq.l7) 


D.  E  Xi  -  a' 

i-Zi 


(Eq.18) 


D.  (b+1)  E 

i-Zi 


(Eq.l9) 


Using  a",  we  can  apply  Che  antiderivative  of  the  power  law  to  integrate  over 
the  unmeasured  regions  of  Che  profile  at  the  cop  and  bottom  (see  1[F-1.6  and 
•[F-1.7). 
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F-1.3.  DETERMINING  NEAR-SHORE  DISCHARGE.  You  can  use  the  TRANSECT  PLAYBACK 
menu  Co  estimate  the  discharge  near  the  shore.  TRANSECT  uses  a  ratio -Interpo¬ 
lation  method  for  estimating  the  velocity  between  the  channel  bank  and  the 
first  or  last  known  mean  velocity  (V^)  in  a  subsection  (ensemble  or  segment  if 
averaging  is  ON  in  PLAYBACK) .  If  we  assume  a  triangular  area  between  the  last 
subsection  and  the  channel  bank,  an  algorithm  for  estimating  the  discharge  in 
this  near  shore  area  is : 


Q.id.  -  0-707  V„  L  d,  /  2.  (Eq.20) 

where:  d,,  -  Depth  of  Che  first  or  last  subsection. 

L  —  Distance  from  the  shore  to  the  last  subsection. 

Again,  note  that  if  averaging  is  ON  in  £LAYBACK,  V„  is  the  mean  velocity  for 
the  segment  (averaging  interval).  If  averaging  is  OFF,  is  Che  mean  veloci¬ 
ty  for  Che  ensemble.  This  allows  flexibility  if  there  are  bad  ensembles  at 
the  start  or  end  of  Che  transect. 

F-1.4.  HOW  TRANSECT  DETERMINES  THE  SIZE  OF  THE  TOP/BOTTOM/MIDDLE  HATER  LAY¬ 
ERS.  Before  computing  discharge,  TRANSECT  must  determine  the  size  of  the  top, 
bottom,  and  middle  water  layers  (see  Figure  F-2).  See  Table  F-1  for  term 
definitions . 

TRANSECT  uses  Che  center  of  Che  first  depth  cell  to  find  the  thickness  of  Che 
cop  water  layer.  To  compute  the  thickness  of  Che  top  water  layer,  we  start 
with  Che  depth  Co  the  center  of  the  first  depth  cell,  D^op. 

Dtop  “  ^ADCP  +  Dj,  -t-  ((Dp  -t-  Dp  +  D,)  /  2)  (Eq.21) 

For  the  total  water  depth,  we  use  the  average  of  the  beam  depths. 

D.ir*  “  Average  of  four  beam  depths  (Eq.22) 

TRANSECT  computes  the  last  good  bin  depth  from  D.v*  by  determining  the  depth 
of  noise  interference  from  the  acoustic  side  lobes,  adding  Che  depth  of  the 
transducer  face  from  the  surface ,  and  adding  a  thickness  chat  depends  on  the 
transmit  pulse  sequence. 

Dir _ -  (D„,  cos(e)  +  Djujcp)  •  ((Dp  -  Dp)  /  2)  (Eq.23) 

D, _ _ is  the  last  possible  depth  of  a  good  bin.  Dlc  is  the  depth  of  the  low¬ 
est  bln  that  is  above  D, _ _  The  position  of  the  last  good  depth  cell  (Dy-) 

gives  us  Che  starting  depth  of  the  bottom  water  layer.  The  valid  ADCF  veloci¬ 
ty  data  in  depth  cells  starting  at  Dppp  and  ending  at  Di,c  are  used  Co  calculate 

Che  middle  layer  discharge  (MidQ) .  TRANSECT  obtains  the  distance  Co  Che  bot¬ 

tom  (Dppppi)  as; 


Dtot.»i  "  Dpvi  +  D^cp  (Eq.24) 

From  the  preceding  we  can  define  the  boundaries  of  the  water  layer  (Fig.  F-2): 

Zi  ”  Dpptpi  -  Di,(;  -  Dp/2  (Eq.25) 

Zj  -  Dppp.i  -  Dppp  -H  D./2  (Eq.26) 

Z,  -  D,pp.,  (Eq.27) 
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The  water  layer  thickness  follows ; 


Top  layer  -  Z3 
Middle  layer  -  Z2 
Bottom  layer  -  Zj 


(Eq.28) 

(Eq.29) 

(Eq.30) 


F-1.5.  HOW  TRANSECT  CALCULATES  HIDDLE  UYER  DISCHARGE  (MidQ).  TRANSECT  cal¬ 
culates  the  middle  layer  discharge  over  the  range  determined  by  the  middle 
water  layer  thickness  (Eq.29).  For  each  bin,  TRANSECT  computes  the  discharge 
in  that  bin  using  the  cross-product  (Eq.6)  and  the  time  difference  between 
successive  ensembles. 


Qi  “  At  Dj 

where:  i  -  Bin  number. 

Qi  -  Discharge  in  the  i'*'  bin. 


(Eq.31) 


The  value  of  is  that  value  displayed  in  the  TABULAR  submenu  at  its  corre¬ 
sponding  depth.  TRANSECT  determines  the  discharge  for  the  middle  layer  by 
summing  all  the  discharges  from  the  individual  bins . 


MidQ  -  Z  Qi  (Eq.32) 

i-1 

where:  i  -  The  first  good  bin.  If  the  first  bin  has  bad  data,  TRANSECT  in¬ 
crements  i  until  it  finds  good  data.  If  TRANSECT  does  not  find 
good  data,  it  does  not  calculate  discharge  for  the  ensemble. 

N  -  The  number  of  the  bln  just  above  the  bin  that  Dls  intersects. 

F-1.6.  HOW  TRANSECT  ESTIMATES  TOP  LAYER  DISCHARGE.  There  are  two  ways  to 
compute  the  top  layer  discharge  estimate  based  on  the  setting  of  TOP  DISCHARGE 
ESTIMATE  in  the  configuration  file.  The  two  options  are  CONSTANT  and  POWER.  Para¬ 
graph  ^F-1.2  outlines  these  methods. 

Constant  Method,  transect  estimates  the  discharge  for  the  top  water  layer  by 
extrapolating  the  value  of  the  cross-product  in  the  first  good  bin  to  the 
surface.  This  method  extrapolates  data  in  a  straight  line  to  the  surface. 


TopQ  "  XpQ  At  (Z3  -  Z2) 


(Eq.33) 


Power  Method,  transect  estimates  the  discharge  for  the  top  water  layer  by 
integrating  the  power  law  of  Eq.  10  over  the  top  water  layer  and  multiplying 
by  At  (see  Figure  F-2): 

Zj 

TopQ  -  At  a"  J  z’>  dz  (Eq.34) 

Z2 


At  D.  (Zj‘^‘  -  Zj"*^)  E 

i-Zi 

Zb+l  7  b'9-1 
2 


(Eq.35) 
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F-1.7.  HOW  TRANSECT  ESTIMATES  BOHOM  LAYER  DISCHARGE.  There  are  two  ways  to 
compute  the  bottom  layer  discharge  estimate  based  on.  the  setting  of  BOTTOM  DIS¬ 
CHARGE  ESTIMATE  In  the  configuration  file.  The  two  options  are  CONSTANT  and  POWER. 
Paragraph  ^F-1.2  outlines  these  methods. 

Constant  Method,  transect,  estimates  the  discharge  for  the  bottom  water  layer 
by  extrapolating  the  value  of  the  cross-product  In  the  last  good  bln  (X^c)  to 
the  bottom.  This  method  extrapolates  data  In  a  straight  line  to  the  bottom. 

BtmQ  -  Xu:  At  Zj  (Eq.36) 

Power  Method.  TRANSECT  estimates  the  discharge  for  the  bottom  water  layer  by 
Integrating  the  power  law  of  Eq.  10  over  the  bottom  water  layer  and  multiply¬ 
ing  by  At  (see  Figure  F-2) : 


Zi 

BtmQ  -  At  a"  j"  z*"  dz 
0 

Zr 

At  D,  Zi"*^  E  Xi 
1-Z: 

Zj!-*!  -  Zi*^^ 


(Eq.37) 


(Eq.38) 
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